Mon. Not. R. Astron. Soc. 000, 000-000 (199?) Printed 7 February 2008 (MN IM^X style file vl.4) 



Accuracy of the relativistic Cowling approximation 
in slowly rotating stars 



On 
On 



On 



oo 

o 
in 
o 

On ■ 
"o ! 

cr 

i ' 



Shijun Yoshida 1 and Yasufumi Kojima 2 

1 Department of Earth Science and Astronomy, Graduate School of Arts and Sciences, University of Tokyo, 
Komaba, Meguro, Tokyo 153, Japan 

2 Department of Physics, Hiroshima University, 
Higashi- Hiroshima 739, Japan 



Accepted ? Received ?; in original form 



ABSTRACT 

We have calculated the non-radial oscillation in slowly rotating relativistic stars with 
the Cowling approximation. The frequencies are compared with those based on the 
complete linearized equations of general relativity. It is found that the results with 
the approximation differ by less than about 20% for typical relativistic stellar models. 
The approximation is more accurate for higher-order modes as in the Newtonian case. 
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1 INTRODUCTION 

In recent years, our understanding of pulsations of non- 
rotating relativistic stars has been much improved. In par- 
ticular, the role of gravitational wave in the stars becomes 
much clear. There exists the oscillation mode named as w- 
mode (gravitational uiave mode), associated with the gravi- 
tational wave. The gravitational wave is inherent in the gen- 
eral relativity, so that the mode becomes evident only for 
relativistic system. Except for this new mode, the relativity 
little affects the modes known in the Newtonian pulsation 
theory. The general relativity slightly changes the oscilla- 
tion frequency and gives rise to a very slow damping of the 
mode. The emission of gravitational radiation implies that 
the oscillation frequency should be complex with a relatively 
t iny im aginary part. (See Andersson, Kojima, and Kokkotas 
(1996) and references therein for the present status of the 
oscillation spectra of non-rotating relativistic stars.) 

The perturbation of the gravity is not so important, 
if the w-mode and the decay of the pulsations due to the 
gravitational radiation do not matter. This fact hints at the 
further simplification of the pulsation problem, that is, ne- 
glecting the perturbation of the gravitational field. This ap- 
proximation is known as the Cowling approximation (Cowl- 
ing (1941)) in the Newtonian stellar pulsation theory, and 
gives the same qualitative results and reasona ble a ccuracy 
of the oscillation frequencies. (See, e.g., Cox (1980).) Two 



different prescriptions, so far, had been proposed for the rel- 
ativistic Cowling approximation in non-radial pulsations of 
the non-rotating stars. One method is th at all metric per- 
turbations are neglec ted, i.e., Sg^ = ( McDermott, Van 
Horn fc Scholl 1983 ). The other is that the Sg r t compo- 
nent of the metric perturbations is retained in the pulsation 



equations ( Finn 1985 ) . Lindblom and Splinter ( 1990] ) exam- 
ined the accuracy of these two versions of the approximation 
for the dipole p-modes, and concluded that the McDermott, 
Van Horn, and Scholl version is more accurate than Finn's 
version. How ever, the F inn's version is superior in g-mode 
calculations ( Finn 198^ ). 

In contrast to the oscillations in the non-rotating stars, 
the calculation of oscillation frequencies in rapidly rotating 
relativistic stars seems to be very difficult task. The equa- 
tions become significant ly com plicated due to the relativity 
and rotation. See Priou (1992) for a set of the explicitly writ- 
ten, lengthy equations. Therefore, the calculation of normal 
frequencies was limit ed to the non -rotating and slowly ro- 
tating stars at best (Kojima 1992 ) . The possibility to use 
the Cowling approximation in rotating r elativ istic system 
should be studied. Ipser and Lindblom (1992) formulated 



the Cowling approximation of the relativistic system, set- 
ting the metric perturbation to zero. If the approximation is 
applied to the stellar pulsation, the equation becomes single 
second order partial differential equation. The equation is 
manageable and will be hopefully solved in the future. The 
crucial point is how good the Cowling approximation is to 
estimate the oscillation frequency in the relativistic regime. 

In this paper, we will estimate the accuracy of the Cowl- 
ing approximation in the frequencies of the non-rotating star 
and its rotational corrections. That is, we will compare the 
frequencies of the Cowling approximation with those of the 
relativistic perturbation theory. In section 2, we present the 
formalism to calculate normal frequencies and the rotational 
shifts by the the Cowling approximation. In section 3, the 
comparison between both numerical results is given. Finally, 
section 4 is devoted to the discussion. Throughout this paper 
we will use units of G — c — 1. 
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2 METHOD 

To assess the accuracy of the Cowling approximation, two 
distinct calculations are performed. In one calculation, the 
relativistic perturbation equations are solved to find nor- 
mal frequencies and the rotational corrections for the pul- 
sation of stellar models. The equations and techniques are 
described in Kojima (1992,1993). In the other calculation, 
we adopt the relativistic Cowling approximation. We restrict 
ourselves to the slowly rotating case, and write down the 
equations of the pulsation and rotational corrections below. 
Normal frequencies are real numbers, because gravitational 
perturbations are neglected in the Cowling approximation. 
On the other hand, the normal frequencies are complex num- 
bers with small imaginary part representing the damping 
of the gravitational radiation, when solving the exact pul- 
sation equations. We make a comparison in the oscillation 
frequency only. 



2.1 Equilibrium state 

We assume that the star is slowly rotating with a uniform 
angular velocity Q. In this paper we consider only the first 
order effect with e = f2/-v/ M/R 3 , where R and M are the 
radius and the mass of the star, respectively. The geometry 
of a slowly rotating star in general relativity is described by 



ds 



+ e x dr 2 + r 2 



e u dt' 



-2ur 2 sin 2 6 dtdcj) 



(de 2 



■ sin 2 6 d(j) 2 ) 



(1) 



where u, A and u are functions of r only. Given an equa- 
tion of state p — p(p), the construction of the slowly rotat- 
ing equi libriu m model is well-known. (For detai ls see, e.g., 
Hartle (Il967|); Chandrasekhar and Miller (11974).) We will 



assume a polytropic equation of state, 



(2) 



where n and k are constants. The function u) is of order e, 
and the others, p, A, p, and p are the same as in the non- 
rotating case. 

2.2 The Cowling approximation 

In this subsection, we will explicitly write down the pulsation 
equation for the slowly rotating stars up to the first o rder o f 
e, using the equations derived by Ipser and Lindblom ( 1992 ). 
We consider only the adiabatic perturbation for the perfect 
fluid, and use the simplest model for the adiabatic index F 
based on the structure of the equilibrium configuration, 



r 



p + p dp 



(3) 



p dp 

The frequencies of the all g-modes are therefore zero in 
this case. The perturbations are given in terms of e _lCTt , 
radial function SV(r), and spherical harmonics Yi m (0, cf>), 
where the azimuthal dependence of the spherical harmonics 
is e lm< ^ in our convention. The oscillation frequency a is the 
value measured in the inertial frame. Co rrespo nding to equa- 
tions (16)-(18) of Ipser and Lindblom (1992), the Eulerian 
changes in the pressure, the density, and the four-velocity 
are written up to the first order of e as 

ae- v/2 (p + p)5VY lm e-^ t , (4) 
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where 

VJ = f2 — LO. 



(9) 



(10) 



The terms with zu are rotational corrections to the oscilla- 
tion in the non-rotating star. 

Since the background geometry is spherically symmet- 
ric, the angular part can be separated into each I, m mode. 
The basic equation to be solved becomes second order ordi- 
nary differentia l equa tion for SV. (See equation (38) of Ipser 
and Lindblom Jl992|).) We have 



d . d 2 ^ 

dr dr 



SV = m (2a SIC - S)5V , 



(11) 



where 

A=(p + p)r 2 e { "- x)/2 , 
d 



B = e 



-v/2 



dr 
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r — e 

dr 
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-J(J+1) (p+p)e (l/+A)/2 ,(13) 
(14) 



S = (p + p)r-V"+ A >/ 2 

e -(3,+A)/2 d 



2ovje~ v + 
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(p + p) e K " — (-cue r I 
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(15) 



where S and are of the order e. 

We will solve the eigen-value problem, equation (|ll| ) by 
standard perturbation method, that is, expand the potential 
SV and the frequency a with respect to the small parameter 
e as, 

SV = SV + 8Vie + -- - , (16) 

a = (To + <xi e H . (17) 

The equation of order e° is for the non-rotating case, 
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— A— + B + a 2 C 
ar ar 



SV = 



(18) 



and the equation of order e is 

SVx + 2<T cri C SVo 



dr dr 



M 

TWr* 



(2a flC- So)5V 



(19) 



where So is defined by equation ( |l5| ) with a — ao- Multi- 
plying SVq to equation (|l9|), and integrating over the star, 
we have the first order rotational correction of the normal 
frequency o\ as 



o\ — rn 



M 



f R S \6V \ 2 dr 
2ao f*C\SV \ a dr 



(20) 



Similar expression can be found in the Newtonian pulsation 
equation, but is written by the Lagrangian displacement, not 
by the Eulerian change of the potential. Note also that the 
correction ( |2o| ) can be determined from the eigen-function 
and eigen-value in the non-rotating case, i.e., SVo and ao. 
Therefore, it is sufficient to calculate SVo an d ao- 

The frequency in the Cowling approximation should be 
compared with the real part of eigen- frequency calculated by 
the exact pulsation equation. The frequency a of the slowly 
rotating star can be parameterized as 



a = or + a R me , 
where 

<Jr — 0"o , a R = ffi/r 



(21) 



(22) 



The subscript R means the real part of the exact normal 
frequency. We will compare results of two calculations in or 
and <Jr. 



2.3 Numerical method 

As mentioned in the last subsection, we need eigen-frequency 
and eigen-function of the pulsations in the non-rotating 
stars. The method to solve equation ([iS]) is well-known and 
straightforward. We will briefly summarize the boundary 
conditions and the numerical method. 

Since the basic equation (|l8| ) is a second order differen- 
tial equation, two boundary conditions should be given for 
SVo a t the center and the surface of the star. The boundary 
condition at the center is that the function SVo should be 
regular near r = 0. This condition can be written by power 
series as 



SVo = r l (v + v 2 r' 2 + ■ 



)• 



(23) 



At the surface, the Lagrangian perturbation of pressure 
should be vanished, Ap = . This condition can be written 
explicitly as 



dSVp 
dr 



2 \- v , . , dr 
o"o e (p + p) -T- 
ap 



1 



dp 
p + p dr 



SV 



(24) 



We can numerically integrate the equation (y_8|) from the 
center with the condition (|23J), and from the surface with 
( [ii] ) to the appropriate interior point. The eigen-value ao 
can be obtained as the consistent frequency so as to match 



functions smoothly at the point. We also numerically calcu- 
late the rotational correction of normal frequency o\ , by ao 
and SVo through the equation (pOj). 

In order to check our numerical code for the relativistic 
Cowling approximation, we calculated the fundamental nor- 
mal mode frequencies having harmonic index / in the range 
2 < / < 5 for t he no n-rotating stellar model used by Ipser 
and Lindblom (1992). Our results were in good agreement 



with their normal frequencies and the difference was of the 
order less than 0.05%. 



3 NUMERICAL RESULTS 

We have calculated normal frequencies for a range of the stel- 
lar models with a polytropic equation of state (Pj) . We adapt 
two different polytropic indices (n — 1.0, 1.5). The compact- 
ness of the star ranges from M/R = 0.05 to M/R = 0.20. 
The oscillation frequencies and the rotational corrections are 
calculated for the f-, pi-, and p2-modes having harmonics in- 
dex I in the range 2 < / < 5. The results for typical stellar 
models are shown in Tables 1-4. Normal frequencies of the 
non-rotating stars and the rotational corrections are com- 
pared. Since there is not so much difference in the polytropic 
index, we hereafter concentrate in the model with n = 1. 

We show the normal frequencies and rotational correc- 
tions of the f-modes as a function of M/R in Fig.l. The 
eigen-frequencies of the non-rotating star, or are given in 
Fig. 1(a), and the rotational corrections a' R are in Fig. 1(b). 
The results with the exact relativistic calculation are rep- 
resented by the dashed line, and those with the Cowling 
approximation are by the solid lines. Figure 1(a) shows that 
two results of the f-mode approach with the increase of M/R, 
that is, the error monotonically decreases. For example, the 
relative error for I — 2 f-mode in M/R = 0.05 model is 
about 30% in the magnitude, but it decreases to 15% in 
M/R = 0.2 model. The relativistic Cowling approximation 
therefore gives better result in the f-mode calculation, as 
equilibrium model becomes more relativistic. Figure 1(b) 
shows that both results are always in good agreement on 
the rotational corrections. The relative error is always less 
than 3%. 

Next, we show results of the pi- and p2-modes in Figs. 2 
and 3, respectively. Compared with the exact calculation, 
the dependence of M/R is similar in the oscillation frequen- 
cies and the rotational corrections. The relative error is al- 
most constant in or, but increases in <j'r with M/R. The 
approximation for these p-modes becomes worse for the rel- 
ativistic stars. However, the discrepancy is always within 
10% error. 

From Figs. 1-3, it is quite clear that the accuracy of 
the Cowling approximation increases with the spherical har- 
monic index, I. This feature is well-known in the Newtonian 



Cowling approximation. (See, e.g. Cox (1980).) The neglect 



of the gravitational perturbation is justified due to the av- 
eraging effect of the local fluctuations for high-order modes. 
The Newtonian Cowling approximation also becomes bet- 
ter wi th the increase of the radial node number. (See, e.g., 



Robe ( 196£ ) for numerical results of the Newtonian Cowling 
approximation.) In the relativistic case, the accuracy also in- 
creases, but the behavior is a little complicated due to the 
additional parameter M/R. 
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4 DISCUSSION 

In this paper, we have compared the relativistic Cowling 
approximation with the perturbation calculation of general 
relativity. The accuracy of the Cowling approximation is 
considerably good, concerning the normal frequencies in the 
non-rotating star and the first order rotational corrections. 
Remarkable fact is that the relative error of the f-mode de- 
creases with the relativistic factor, M / R. The approximation 
is also better for larger radial node number and harmonic in- 
dex I, and greater central condensation, as in the Newtonian 
pulsation theory. 

One of the important applications of the normal modes 
is to determine the critical angular velocity of the rotating 
stars. The maximum angular velocity is limited by the ex- 
istence of the axisymmetric equilibrium state. The equilib- 
rium state however suffers secular instability due to viscosity 
and/or gravitational radiation reaction (GRR). The insta- 
bility sets in through a = mfi (viscous instability) or a = 
( GRR instability), respectively (e.g., Friedman and Schutz 
(1978)). Form the Newtonian calculation, most significant 
contributions to the GRR instability are the f -mode s with 
/ = — m ~ 4. (See, e.g., Yoshida and Eriguchi (1995).) The 
present study shows that the relativistic Cowling approxi- 
mation is able to give a good prediction of the oscillation fre- 
quency within 6% error for the f-modes with I — —m = 4, 5 
in the typical neutron stars, M/R ~ 0.2. The critical point 
is however not within the regime of the linear extrapolation 
from the non-rotating case, as suggested in the Newtonian 
calculation. That is, the present results, an and a' R are not 
sufficient, and the higher order rotational corrections are 
necessary. Therefore, solving normal frequency in rapidly 
rotating stars with the Cowling approximation would be a 
good method to determine the critical angular velocity. 



Table 1. The normal frequencies of the non-rotating stars and 
the rotational corrections for n = 1.0 and M/R = 0.1. Subscripts 
"ex" and "Cow" denote the results with the exact relativistic 
perturbation calculation and the Cowling approximation, respec- 
tively. All entries are given in units of (M/R?) 1 / 2 . 



1 


mode 


C" ex 


& Cow 




a' 

O OW 


2 


f 


1.21 


1.52 


0.577 


0.589 




Pr 


3.11 


3.43 


0.880 


0.883 




P2 


4.oo 


5.09 


0.925 


0.925 


3 


f 


1.62 


1.82 


0.720 


0.723 




Pi 


3.61 


3.85 


0.906 


0.908 




P2 


X. Oft 


5.58 


0.935 


0.936 


4 


f 


l.yo 


2.07 


0.789 


0.791 




Pi 


44. UO 


4.22 


0.923 


0.924 




P2 


O.o4 


6.01 


0.944 


0.944 


5 


f 


Z. ±0 


2.29 


0.831 


0.832 




Pi 


A AC\ 


4.55 


0.934 


0.935 




P2 




6.42 


0.950 


0.951 


Table 2. The 


same 


as Table 1, but for n = 1. 


I 


mode 


C" ex 


&Cow 




a' 

O OW 


2 


f 


1.17 


1.35 


0.670 


0.680 




Pi 


2.70 


2.99 


0.870 


0.878 




P2 


4.12 


4.38 


0.898 


0.900 


3 


f 


1.52 


1.64 


0.780 


0.782 




Pi 


3.19 


3.39 


0.895 


0.900 




P2 


4.64 


4.84 


0.909 


0.911 


1 


f 


1.78 


1.87 


0.832 


0.834 




Pi 


3.59 


3.74 


0.912 


0.915 




P2 


5.10 


5.26 


0.919 


0.920 


5 


f 


2.01 


2.08 


0.865 


0.865 




PI 


3.94 


4.06 


0.925 


0.927 




P2 


5.52 


5.65 


0.927 


0.928 



1.0 and M/R = 0.2. 
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Table 3. The same as Table 1, but for n = 1.5 and M/R = 0.1. 



1 


mode 


<j ex 


&COW 


a' ex 


Cow 


2 


f 


1.43 


1.68 


0.591 


0.609 




Pi 


2.87 


3.18 


0.851 


0.855 




P2 


4.27 


4.56 


0.911 


0.911 


Q 
O 


f 
I 


l.oo 


1 Q7 

i.y * 


U. ( zy 


n *7Q/i 

U. I o4 




PI 


3.35 


3.55 


0.886 


0.889 




P2 


4.76 


4.98 


0.926 


0.926 


4 


f 


2.11 


2.21 


0.796 


0.798 




Pi 


3.73 


3.88 


0.908 


0.910 




P2 


5.19 


5.36 


0.937 


0.938 


5 


f 


2.34 


2.41 


0.835 


0.837 




Pi 


4.Uo 


4.16 


0.923 


0.925 




P2 


5.56 


5.69 


0.945 


0.946 


Table 4. The 


same 


as Table 1, but for n = 1. 


I 


mode 


u ex 


""COW 


ex 


a' 

O OW 


2 


f 


1.36 


1.47 


0.698 


0.706 




Pi 


2.46 


2.71 


0.855 


0.867 




P2 


3.54 


3.82 


0.895 


0.900 


3 


f 


l.UO 


1.75 


0.795 


0.797 




PI 


2.93 


3.08 


0.889 


0.894 




P2 


4.05 


4.24 


0.909 


0.912 


4 


f 


1.93 


1.98 


0.843 


0.844 




Pi 


3.29 


3.40 


0.910 


0.912 




P2 


4.47 


4.61 


0.921 


0.923 


5 


f 


2.13 


2.17 


0.872 


0.873 




Pi 


3.60 


3.68 


0.924 


0.925 




P2 


4.84 


4.95 


0.930 


0.932 



5 and M/R = 0.2. 
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Figure caption 



Fig. 1(a). Normal frequencies, <tr vs. relativistic factor M/R 
for f-mode and n = 1.0. Dashed lines and solid lines corre- 
spond to results with the relativistic perturbation calcula- 
tion and the Cowling approximation, respectively. Normal 
frequencies are normalized by (M/R 3 ) 1 / 2 . Attached labels 
I = 2, 3, 4 and 5 mean the harmonic indices. 



Fig. 1(b). Rotational corrections, a' R vs. relativistic factor 
M/R for f-mode and n = 1.0. Dashed lines and solid lines 
correspond to results with the full calculation and the Cowl- 
ing approximation, respectively, rotational corrections are 
normalized by (M/R 3 ) 11 ' 2 . Attached labels I = 2, 3, 4 and 5 
mean the harmonic indices. 



Fig. 2(a). The same as Fig. 1(a), but for pi-mode 



Fig. 2(b). The same as Fig. 1(b), but for pi-mode. 



Fig. 3(a). The same as Fig. 2(a), but for p2-mode. 



Fig. 3(b). The same as Fig. 2(b), but for p2-mode. 
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Fig.l(b) 




.05 .1 .15 .2 .25 

M/R 



Fig.2(a) 




.05 .1 .15 .2 



M/R 



Fig.2(b) 




.05 .1 .15 .2 .25 



M/R 



Fig.3(a) 




.05 .1 .15 .2 

M/R 



Fig.3(b) 




.05 .1 .15 .2 .25 



M/R 



